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Using lattice simulations of quenched QCD we estimate the finite size effects present when a gluon 
plasma equilibrates in a slab geometry, i.e., finite width but large transverse dimensions. Significant 
differences are observed in the free energy density for the slab when compared with bulk behavior. 
A small shift in the critical temperature is also seen. The free energy required to liberate heavy 
quarks relative to bulk is measured using Polyakov loops; the additional free energy required is on 
the order of 30 - 40 MeV at 2 - 3 T c . 
PACS numbers: ll.10.Wx 12.38. Gc 12.38. Mh 



I. INTRODUCTION 

The formation of a quark-gluon plasma in a central 
heavy ion collision is generally assumed to take place in 
a coin-shaped region roughly 1 fcrmi in width, with ra- 
dius comparable to the radii of the colliding nuclei, which 
is to say several fermi. While lattice gauge theory has 
given us information about bulk thermodynamic behav- 
ior, finite size effects have up to now been studied using 
simplified, phenomenological models. In this work, we 
study via lattice gauge theory simulations the behavior 
of a gluon plasma restricted to a slab geometry, with 
the longitudinal width much smaller than the transverse 
directions. This inner region is heated to temperatures 
above the bulk deconfincmcnt temperature, surrounded 
by an outer region which is kept at a temperature below 
the deconfinement temperature. This outer region pro- 
vides confining boundary conditions for the inner region. 
From our simulations, we derive information about the 
thermodynamics of the gluon plasma in a slab geometry, 
the effective non-equilibrium surface tension and the free 
energy cost of creating quarks relative to bulk. 

Measurements of the equilibrium surface tension cxq of 
pure SU(3) lattice gauge theory (quenched lattice QCD) 
show that the dimensionless ratio cyq/T^ is small. For 
the case of an iV 5 3 x 4 lattice, a /T^ « 0.0292(22). § A 
simple estimate of surface tension effects on the transition 
can be obtained from a simplified model in which only 
volume and surface terms appear, as in the bag model. 0] 
Since we are interested in the pure gauge theory, we take 
the free energy density of the hadronic phase to be zero, 
neglecting the contribution of glueballs, whose masses are 
large compared to the deconfinement temperature. The 
free energy of a gluon plasma of volume V and surface 
area A we take to be given by 



F= f B-18— T 4 ) V mu.1 



(1) 



The bag constant B, which here is simply related to the 
latent heat at the deconfinement transition, is taken to 
be B = (200MeV) 4 . Taking a /T 3 as above, the bulk 
critical temperature to be T c = 260MeV r and the slab 
width as w = Ifm. leads to an increase in T c of about 
3%. Naively, finite size effects due to the surface tension 
should be small. 

Other contributions to finite size effects come from 
a variety of sources. In the case of systems with non- 
abelian symmetries, global color invariance produces an 
additional finite volume effect which will not be consid- 
ered here. j| [|J In general, finite size effects lead 
to a rounding of the transition. |(| This can be taken 
into account in the bag model by a Maxwell construc- 
tion, leading to mixed phases and a broadened critical 
region. A recent treatment for the quark gluon plasma 
can be found in H . We have attempted to avoid these fi- 
nite volume effects by making the transverse dimensions 
large. 



II. METHODOLOGY 

In lattice calculations, finite temperature is introduced 
by the choice of N t , the extent of the lattice in the 
(Euclidean) temporal direction. The relation of phys- 
ical temperature T to Nt and the lattice spacing a is 
simply T — \/N t a. The lattice spacing a implicitly de- 
pends on the gauge coupling j3 in a way determined by 
the renormalization group equations. To lowest order in 
perturbation theory, the relation is given by 
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where is renormalization group invariant and the 
renormalization group coefficients &o and b\ are given by 

117V 34 / N \ 2 

In analyzing our data, we used the renormalization 
group results given in reference ||], which are deter- 
mined directly from lattice simulations, and contain 
non-perturbative information about the renormalization 
group flow. 

By allowing the coupling constant j3 = 6/g 2 to vary 
with spatial location, a spatially dependent temperature 
can be introduced. We have chosen the temperature in- 
terface to be sharp, in such a way that the lattice is di- 
vided into two spatial regions, one hotter and one colder. 
Figure 1 shows the slab geometry that was used for sim- 
ulations. By varying the spatial and temporal plaque- 
tte couplings independently, it is possible in principle 
to make the spatial lattice spacing constant in physi- 
cal units, at least in a slab geometry. || However, the 
quenched approximation simplifies the role of the cold 
region. In the quenched approximation below T c , the 
dominant excitation at low energies is the scalar glue- 
ball. Temperatures near T c are smaller than glueball 
masses by about a factor of four, so glueballs play no 
essential role in the thermodynamics, and the pressure 
in the hadronic phase is essentially zero. For this rea- 
son, we expect the slab thermodynamics to be largely 
insensitive to the precise temperature of the region out- 
side the slab, as long as it is sufficiently low. The outer 
region thus merely provides boundary conditions for the 
slab. A numerical check of this assumption is discussed 
below. In full QCD, this insensitivity to the outer tem- 
perature would not hold, due to pions. Note that the 
role of boundary conditions here is quite different from 
those relevant for the QCD phase transition in the early 
universe. Assuming the transition in full QCD is first 
order, nucleation will occur as the universe cools below 
the transition. If the equilibrium surface tension is small, 
the nucleation rate will be high, and no significant un- 
dercooling is expected. In this spherical geometry 
is relevant, and simulations have been performed where 
Pi n and (3 out have been taken to be near T c . [jl0|, |ll]] 
Note that such simulations are not directly relevant to 
heavy ion collisions, since substantial temperature gradi- 
ents will develop by the time the expanding quark-gluon 
plasma becomes approximately spherical. 

Another difference between this system and near- 
equilibrium systems lies in fluctuations of the interface. 
For equilibrium or near-equilibrium systems, interfacial 
fluctuations play a significant role. In three dimensions, 
we expect that an equilibrium planar interface would al- 
ways be in the rough phase. Interfacial fluctuations as- 
sociated with roughening lead to a universal subleading 
correction to the planar interfacial free energy. [pL2| | , |l3| ] 
In the case of nucleation, interfacial fluctuations play an 
important role in the prefactor of the nucleation rate via 



the functional determinant. |14| In the model we are con- 
sidering, the interface is pinned by the sharp temperature 
interface, which may drastically reduce fluctuations. It 
seems likely to us that this is a physical effect. Further- 
more, the surface modes associated with roughening are 
long-wavelength modes which are likely to require sub- 
stantial time to equilibrate. We conclude that interfacial 
fluctuations are unlikely to play the same role in an an 
expanding quark-gluon plasma created by a heavy ion 
collision that they would in an equilibrium situation. 

In our simulations, we work with slabs of fixed lattice 
width, w — 6a, rather than of fixed physical width. Since 
N t — 4, wT is fixed at 3/2. Maintaining a fixed width 
in unphysical lattice units means that the width of the 
slab in physical units varies by almost a factor of 3 over 
the range of j3 values used, from approximately 1.2 fermi 
to 0.5 fermi. At higher temperatures, this is somewhat 
smaller than the longitudinal size of the plasma formation 
region expected in heavy ion collisions. 

The use of equilibrium statistical mechanics to study 
gluon plasma properties during the early stages of plasma 
formation may appear suspect. It is useful to compare 
the conditions here with those of the Bjorken model. fjl| 
In this model, the temperature is specified on an initial 
surface of constant proper time To, and has a simple scal- 
ing behavior given by 

T(t)=T (t /t/s (4) 

w here c s is the speed of sound. For an ideal gas, c s — 
This can be used to estimate the spatial varia- 
tion of the temperature of a quark-gluon plasma shortly 
after formation in a central collision. For example, when 
a coin-shaped region of width 1 fermi has expanded to 
1.5 fermi, the variation in temperature is only from O.8T0 
at the center of the coin to T at its edges. It has been 
argued that c s may be much smaller than \/l/3 near T c 
p6[ , p7[ ], which would further decrease variation in T. 
Thus the assumption of constant temperature is likely 
to be a good approximation to the actual early stages 
of a quark-gluon plasma formed in a heavy ion colli- 
sion. Although, a non-trivial temperature profile could 
be incorporated by adjusting the spatial dependence of 
the coupling constants, it is neither necessary nor desir- 
able. A more significant difference between the condi- 
tions studied here and those of the Bjorken model lies in 
the velocity distribution. The Bjorken model leads to a 
space-time dependent velocity distribution whereas lat- 
tice models are restricted to an average velocity of zero. 
Adding even a constant average velocity has a technical 
difficulty similar to the well-known problem of adding a 
non-zero chemical potential to lattice models. 

The free energy density / for the slab was obtained 
using the standard method Q of integrating the lattice 
action with respect to (5. We use a convenient convention 
for the sign of / that is opposite the usual one. In the 
bulk case, / is then identical to the pressure p. 
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f4lL t =Ar * 4 / df3'[(S) T -(S) ] (5) 

where (S) is the expectation value of a plaquette averaged 
over the region of interest: 

S = ^ReTr U p (6) 

The subscripts T and denote expectation values mea- 
sured at finite and zero temperature, respectively. As 
in the bulk case, it is necessary to subtract the zero- 
temperature action density expectation value from the 
finite temperature expectation value, removing terms 
which would give divergences in the continuum limit. 
The subtracted term is obtained from a zero tempera- 
ture simulation with the same pair of (3 values. This is 
clearly necessary because the introduction of two regions 
with differing (3 values introduces an interface even at 
zero temperature. However, there is the possibility that 
new divergences are introduced by the interface between 
the two regions. Such divergences would manifest in the 
need for counterterms on the boundary between the two 
regions. In general, quantum field theories with bound- 
aries develop divergences that are not present in infinite 
volume or with periodic boundary conditions. Symanzik 
has shown to all orders in perturbation theory that 
in the case of <j) A with so-called Schrodinger functional 
boundary conditions, the theory is finite in perturbation 
theory after adding all possible boundary counterterms 
of dimension d < 3 consistent with the symmetries of the 
theory. It is generally believed that this result applies 
as well to all renormalizable field theories and general 
boundary conditions, but a proof is lacking. Luscher et 
al. [^0| have shown for gauge theories that at one loop 
no new divergences are introduced by Schrodinger func- 
tional boundary conditions. This is consistent with the 
non-existence of gauge-invariant local fields of dimension 
< 3 in pure Yang-Mills theory. 

In order to take advantage of the data on bulk thermo- 
dynamics provided by the Bielefeld group |8j , we worked 
consistently with lattices of overall size IfFx 4. Up to 
8000 initialization sweeps were used, and up to 30,000 
sweeps were used for measurements. Measurements were 
made every 10 sweeps, and corrected for autocorrela- 
tion. The values used for each subtraction come from 
16 4 lattices with identical values of f3i n and (3 ou t- The 
value of (3 ou t was held fixed at 5.6 while (3i n varied 
from 5.6 to 6.3. For comparison, the bulk transition 
for N t = 4 occurs at c (N t = 4, N s = 16) = 5.6908(2) 
l3 c (N t = 4, N s = oo) = 5.6925(5). | 

III. FREE ENERGY OF GLUONS 

Figure 2 shows the free energy density //T 4 versus 
T/T c compared with the bulk pressure. The free en- 
ergy in the slab is lower than the bulk value by almost 



a factor of two at 2T C . It appears that the slab value is 
slowly approaching the bulk value, but other behaviors 
are also possible. Calculations of the finite-temperature 
contribution to the Casimir effect for a free Bose field 
contained between two plates show that //T 4 has a non- 
trivial dependence on the dimensionless combination wT . 
pl[ , |2^] It is natural to ask if the corrections to the free 
energy seen here can be accounted for by the conven- 
tional Casimir effect. A straigtforward calculation of the 
free energy of a non-interacting gluon gas confined to a 
slab shows an increase in the free energy density over 
the bulk value by a factor of about 1.63 at wT = 3/2. 
The Casimir effect alone is thus unable to explain the 
reduction of the free energy observed in our simulations. 
We are currently investigating more elaborate theoretical 
models which include the effect of a non-trivial Polyakov 
loop. 

A consistency check was performed on the surface ef- 
fects using a method originally developed for measuring 
the equilibrium surface tension. ]23( ] The free energy den- 
sity was calculated for a system at /3j„ = 6.0 by perform- 
ing simulations with (3i n fixed at 6.0 and (3 ut varying 
from 5.6 to 6.0. Combining these results with the bulk 
data of reference Q creates a path equivalent to vary- 
ing ftin while holding /3 out fixed. Figure 3 shows the 
two equivalent paths. For /3j„ = 6.0 and (3 ou t = 5.6, 
this gives //T 4 = 0.65 ± 0.04, to be compared with 
f/T 4 = 0.69 ± 0.03 for the direct calculation. 

The major source of systematic error lies with the 
choice of boundary conditions for the slab, here set by 
Pout- We have estimated the effects of varying (3 ou t by 
performing simulations at /3,-„ = 6.2 and f3 out = 5.5 on 
16 3 x 4 and 16 4 lattices. These results suggest that lower- 
ing (3 out from 5.6 to 5.5 reduces the free energy by roughly 
10 percent at /?;„ = 6.2. 

IV. SURFACE TENSION 

The equilibrium surface tension ao can be defined 
from the excess free energy as (3i n and (3 ou t approach 
/3 C , the bulk critical value, with the width w simultane- 
ously taken to infinity. It is common to assume that this 
quantity can also be used to characterize non-equilibrium 
conditions as well. As we have seen in the introduction, 
the small equilibrium value of ao /T 3 would lead to neg- 
ligible surface effects in QCD. 

We define an effective surface tension a(w, T) by 

f=p-2a(w,T)/w (7) 

where the notation a(w,T) recognizes that the surface 
tension a does depend on the width of the slab and the 
internal and external temperatures. The factor of 2 oc- 
curs because the slab has two faces. In the limits where 
T approaches T c and w goes to infinity, this quantity ap- 
proaches ao. Figure 4 shows a(w, T)/T 3 versus T/T c for 
wT = 3/2; representative error bars are shown. 
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The value at j3 = 5.7, 0.056 ± 0.002 is substantially 
higher than the value a = 0.0292 ± 0.0022 given in 
reference We attribute this to two effects: in our 
case p ou t is fixed at 5.6, whereas for equilibrium measure- 
ments it is extrapolated to f3 c , and our finite value of the 
width w also acts to increase a(w, T) over the equilibrium 
value as measured in simulations at large w. Away from 
the bulk critical point, a/T 3 rises quickly to a peak at 
about 1.4T C , and then falls slowly as T increases. A large 
non-equilibrium surface tension has also been observed in 
measurements of the equilibrium surface tension, where 
these effects were obstacles to obtaining ao- p3| 



change occurs between (3 = 5.8 (T = 1.23 T c ) and 
(3 = 5.85 (T = 1.36 T c ). For larger values of f3, AF Q di- 
minishes to a value of approximately 30 — 40 MeV in the 
middle of the slab. In table one, we list (3, T, slab width 
in fermis, width of the core in fermis, and AFq in MeV 
for representative values. The width of the core is cal- 
culated by interpolating the Polyakov loop profiles and 
determining the region where the slab expectation value 
is greater than 80% of the bulk value. All conversions to 
physical units are performed by taking the string tension 
a to be (425 MeV) 2 , which implies T c = 254 MeV. @ 



V. FREE ENERGY OF QUARKS 

The Polyakov loop defined by 

,1/T 

P(x) = (l/N c )Tr Pexp i / A (x,r) elr 

Jo 



(8) 



is the order parameter for the deconfincmcnt transition 
in pure (quenched) gauge theories. In the case of SU(N), 
there is a global Z(N) symmetry which ensures at low 
temperature that the expectation value (TrP) is 0. At 
sufficiently high temperatures, this symmetry is sponta- 
neously broken. The expectation value of the Polyakov 
loop can also be interpreted in terms of the free energy 
of an isolated, infinitely heavy quark Fq: 



{P{x))=ex V [-F Q {S)/T] 



(9) 



In the low-temperature confined phase, Fq is taken to 
be infinite, whereas in the high temperature phase it is 
finite. 

Direct extraction of Fq from computer simulations is 
problematic, because the expectation value has a mul- 
tiplicative, /3-dependent ultraviolet divergence. This di- 
vergence can be eliminated when comparing bulk expec- 
tation values to those in finite geometries. We define 



AF Q (x) = -T ln[P s i ab {x)/Pi 



bulk] 



(10) 



as the excess free energy required to liberate a heavy 
quark in the slab geometry relative to bulk quark matter 
at the same temperature. This technique can also be 
used in, e.g., a spherical geometry, which is relevant for 
nucleation. gj, 

Unlike the thermodynamics, the Polyakov loop pro- 
file is expected to be sensitive to the outer temperature, 
since the Polyakov loop will decay away from the inter- 
face as exp[— err /Tout], where a is the string tension. In 
figure 5, we show the expectation value for the Polyakov 
loop versus z measured in lattice units for several val- 
ues of (3. Each curve is normalized by dividing the val- 
ues of the Polyakov loop by the bulk expectation value 
at the corresponding value of f3. Error bars are shown 
only for even values of z. It is clear that a significant 



VI. CONCLUSIONS 

There are significant deviations in the slab geome- 
try from bulk behavior and the ideal gas law, arising 
from a strong non-equilibrium surface-tension. This non- 
equilibrium surface tension can be an order of magnitude 
greater than the equilibrium value. Surface tension ef- 
fects also produce a mild elevation of the apparent crit- 
ical temperature. Measurement of Polyakov loop expec- 
tation values relative to bulk shows that the suppression 
of heavy quark production due to the slab geometry is 
small. 

There is little doubt that most of the effects seen can 
be attributed to the smallness of the width w relative to 
the natural interfacial width. If w were made arbitrarily 
large, the two interfaces on either side of the slab would 
be independent, and each would have some thickness t, as 
measured, say, from the behavior of the Polyakov loop. 
In this situation, the thin wall approximation, familiar 
in nucleation theory, is valid, and the free energy can be 
decomposed into volume and surface terms. J25J As w is 
reduced and becomes commensurate with t, the thin wall 
approximation breaks down. Intuitively, it becomes im- 
possible to fit two interfaces plus a bulk region between 
them into the width available. For any fixed temper- 
ature and sufficiently small w, no true bulk phase can 
be formed. This is a physical concern for quark-gluon 
plasma formation in heavy ion collisions, independent of 
any computer simulation or model. 

There are good reasons to call our result an estimate 
rather than a calculation. Although lattice gauge theory 
simulations of bulk behavior can be made arbitrarily ac- 
curate in principle, in this case there is some uncertainty 
in the precise exterior boundary conditions appropriate, 
and indeed in the applicability of equilibrium thermody- 
namics at this early stage of quark-gluon plasma forma- 
tion. However, this seems like the best estimate available 
now, and further refinements are possible. An alternative 
approach might be based on real-time simulations; while 
this area is the subject of active research, it is still in its 
infancy. {§, @ 

We have not yet explored the nature of the phase tran- 
sition, which will require some care. One interesting pos- 
sibility is that the order of the transition might change as 
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the width changes. The deconfinement transition in bulk 
quenched finite temperature QCD is in the universality 
class of the three-dimensional three-state Potts model, 
which has a first-order phase transition. As the width 
of the slab becomes commensurate with the correlation 
length near T c , the phase transition should cross over to 
the universality class of the two-dimensional three-state 
Potts model. The two-dimensional three-state Potts 
model has a second-order phase transition, so it is possi- 
ble that the order of the transition may change. |2g] The 
correlation length at the bulk transition is known to be 
large p9j , so it is likely that the transverse correlation 
length in the gluonic sector is much larger in the slab 
geometry than in bulk, even if crossover does not take 
place. 
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FIG. 1. Schematic Drawing of slab geometry. 
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FIG. 2. Free energy density f/T 4 versus T/T c for the slab 
geometry at wT — 3/2 compared with bulk behavior. The 
bulk data is from reference tel. 
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FIG. 3. Paths in the Pin-Pout plane. 
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FIG. 4. Effective surface tension a/T 3 versus T/T c for the 
slab geometry at wT = 3/2. 
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FIG. 5. Polyakov loop expectation value versus position in 
the z direction. 
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